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Using the eigen-mode of the QCD Dirac operator D = y^D^, we develop a manifestly gauge- 
covariant expansion and projection of the QCD operators such as the Wilson loop and the 
Polyakov loop. With this method, we perform a direct analysis of the correlation between con- 
finement and chiral symmetry breaking in lattice QCD Monte Carlo calculations. Even after 
removing the low-lying Dirac modes, which are responsible to chiral symmetry breaking, we 
find that the Wilson loop obeys the area law, and the string tension or the confinement force is 
almost unchanged. We find also that the Polyakov loop remains to be almost zero even without 
the low-lying Dirac modes, which indicates the Z3-unbroken confinement phase. These results 
indicate that one-to-one correspondence does not hold between confinement and chiral symmetry 
breaking in QCD. 
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1. Introduction: relation between confinement and chiral symmetry breaking 

Color confinement and chiral symmetry breaking [|J are striking nonperturbative phenom- 
ena of quantum chromodynamics (QCD). To clarify their correspondence is an important subject 
[Q, |3[ Q|, however, their relation is not yet clarified directly from QCD. The strong correlation be- 
tween them has been suggested by the simultaneous phase transitions of deconfmement and chiral 
restoration in lattice QCD both at finite temperature [Q| and in a small- volume box [Q]. 

The close relation between confinement and chiral symmetry breaking has been also suggested 
in terms of the monopole degrees of freedom [[| [3j, which topo logically appears in QCD by taking 
the maximally Abelian (MA) gauge [§, 0, [|||. Actually, by removing the monopoles, confinement 
and chiral symmetry breaking are simultaneously lost in lattice QCD [J3j] , as schematically shown in 
Fig. 1. This indicates an important role of the monopole to both confinement and chiral symmetry 
breaking, and these two nonperturbative QCD phenomena seem to be related via the monopole. 




Figure 1: The role of monopoles to nonperturbative QCD. In the MA gauge, QCD becomes Abelian-like due 
to the large off-diagonal gluon mass of about lGeV and monopole current topologically appears [0, |[]. 
By the Hodge decomposition, the QCD system can be divided into the monopole part and the photon part. 
The monopole part has confinement [g], chiral symmetry breaking ^ and instantons JlO[], while the photon 
part does not have all of them. Thus, lattice QCD studies indicate the essential contribution of monopoles to 
both confinement and chiral symmetry breaking. However, the direct relation between them is unclear. 



As a possibility, however, to remove the monopoles may be "too fatal" for nonperturbative 
properties. If this is the case, nonperturbative phenomena are simultaneously lost by their removal. 

In fact, if only the relevant ingredient of chiral symmetry breaking is carefully removed, how 
will be quark confinement? To obtain the answer, we perform a direct investigation between con- 



finement and chiral symmetry breaking, using the Dirac-mode expansion and projection [11]. 

2. Gauge-invariant formalism of Dirac-mode expansion and projection 

In this paper, using the eigen-mode of the QCD Dirac operator E> = y^D^, we propose a 
manifestly gauge-covariant expansion/projection of QCD operators such as the Wilson loop and 
the Polyakov loop, and study the relation between confinement and chiral symmetry breaking [O]. 
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2.1 Eigen-mode of Dirac operator in lattice QCD 

In lattice QCD with spacing a, the Dirac operator D = y^D^ is expressed with U^(x) as 



1 4 

®*,y E rE/n [Un(x)$x + d, y - U^(x)8 x ^ y ] , (2.1) 

La jll=l 



with U-n (x) = Ua (x — p.). Adopting hermitian /-matrices = y^ , D is anti-hermitian and satisfies 
Vr y x = —S) x ,y The normalized eigen-state \n) of the Dirac operator D is introduced as 

TD\n) =iXn\n) (2.2) 



with A„ € R. Because of {75,®} = 0, the state 75 \n) is also an eigen-state of E) with the eigenvalue 
— ikn. The Dirac eigenf unction \\f n {x) = (x\n) obeys By n {x) = iX n y n {x), and its explicit form of 
the eigenvalue equation in lattice QCD is 

1 4 

tt L YpPh ( x ) V* ( x + A) - U-n (x) \if„ (x - fL)] = iX n Yn (x). (2.3) 
La \l=\ 

The Dirac eigenfunction \\f n {x) can be numerically obtained in lattice QCD, besides a phase factor. 
According to U^(x) — > V (x)U^(x)V^ (x + jl), the gauge transformation of Yn{x) is found to be 

\l/n(x)^V(x)Yn(x), (2.4) 

which is the same as that of the quark field. To be strict, for the Dirac eigenfunction, there can 
appear an irrelevant ^-dependent global phase factor as e l(p ^ v \ according to the arbitrariness of the 
definition of y n (x). 

From the Banks-Casher relation [12], the quark condensate {qq), the order parameter of chiral 
symmetry breaking, is given by the zero-eigenvalue density p(0) of the Dirac operator D: 

(qq) = - lim lim np(0), (2.5) 

m— >0V— 

where the spectral density p(A) is given by p(A) = yY,n(S(h ~ Ki)) with space-time volume V. 
Thus, the low-lying Dirac modes can be regarded as the essential modes responsible to spontaneous 
chiral-symmetry breaking in QCD. 

2.2 Operator formalism in lattice QCD 

The recent analysis of QCD with the Fourier expansion of the gluon field quantitatively reveals 
that quark confinement originates from low-momentum gluons below about lGeV in both Landau 
and Coulomb gauges [13]. This method seems powerful but accompanies some gauge dependence. 
To keep the gauge symmetry manifestly, we take the "operator formalism" in lattice QCD [11]. 

We define the link-variable operator by the matrix element of 

(x\Ufi\y) = C/ M (x)4+A,y ( 2 - 6 ) 
The Wilson-loop operator W is defined as the product of along a rectangular loop, 

N 

W = l\U l x k = tJ^% 2 ---% N . (2.7) 

k=i 
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For arbitrary loops, one finds Yjk=i Afc = 0- We note that the functional trace of the Wilson-loop 
operator W is proportional to the ordinary vacuum expectation value (W) of the Wilson loop: 

TrW = trj>|#|*) =trY,(x\U lll U lJa ---t) llN \x) 

X X 

= tr (xi\U^\x2}(x 2 \U^ 2 \x 3 )(x 3 \U tl3 \x 4 .)---{x N \U^ N \xi) 

2 N-l 

= tr^{x\U^\x + fa){x + fa\U^\x+Y,fa)---(x + £ A*|*Vl*) 

x k=l k=l 

2 N-l 

= ^{Un 1 (x)U fl2 (x+fa)U fl3 (x+Y,N---U flN (x+Y, fa)} = {W)-Trl. (2.8) 

-v k=l k=l 

Here, "Tr" denotes the functional trace, and "tr" the trace over SU(3) color index. 

The Dirac-mode matrix element of the link- variable operator can be expressed with ljf n (x): 

(m\U\n) = Y,(m\x)(x\U fl \x+fi){x + £i\n) = Y, x l f m(x) u ii{ x )Yn(x+ fa. (2.9) 

X X 

Although the total number of the matrix element is very huge, the matrix element is calculable and 



gauge invariant, apart from an irrelevant phase factor. Using the gauge transformation ( |2.4| ), we 
find the gauge transformation of the matrix element as [pi]] 



{m\Un\n) = Y,Yl( x )U»(x)Yn(x + fa 

X 

-> Z x > / m(x)V t (x)-V(x)U^x)v\x + fa-V(x + fa ¥n (x + fa 

X 

= Y,wi{x)U )1 {x)\if n {x+ fa = (m\0p\n). (2.10) 

x 

To be strict, there appears an ^-dependent global phase factor, corresponding to the arbitrariness of 
the phase in the basis \n). However, this phase factor cancels as e~" Pn e" p " = 1 between \n) and (n|, 
and does not appear for QCD physical quantities including the Wilson loop. 

2.3 Dirac-mode expansion and projection 

From the completeness of the Dirac-mode basis, \n)(n\ = 1, arbitrary operator O can be 
expanded in terms of the Dirac-mode basis \n) as 

d = YY\n)(n\6\m)(m\, (2.11) 



which is the theoretical basis of the Dirac-mode expansion [llj. Note that this procedure is just the 
insertion of unity, and is of course mathematically correct. 

Based on this relation, the Dirac-mode expansion and projection can be defined. We define the 
projection operator P which restricts the Dirac-mode space, 

P=£\n){n\, (2.12) 

ne A 



4 



Confinement and chiral symmetry breaking with Dirac-mode expansion 



Hideo Suganuma 



where A denotes arbitrary set of Dirac modes. In P, the arbitrary phase cancels between \n) and (n\. 

A/} A _r_ A 

One finds P = P and P f = P. The typical projections are IR-cut and UV-cut of the Dirac modes: 

Pm= £ \n){n\, P VY = £ \n}(n\. (2.13) 
W>Air |A„|<A UV 

Using the projection operator P, we define the Dirac-mode projected link-variable operator, 
U p = PU^P= £ £ (2.14) 

During this projection, there appears some nonlocality in general, but it would not be important 
for the argument of large-distance properties such as confinement. From the Wilson-loop operator 
W = IXt=i we define the Dirac-mode projected Wilson-loop operator W = Y\k=i Uu k , and 
rewrite its functional trace in terms of the Dirac basis as [p|] 



N 

r> ■ ■ AD A D A D AD A y\ ,«.,«. ^. A. A, A, 

Tr W p = Tr]J U P k = Tr [/£ £/£••• t/* = Tr PU^PU^P ■ ■ ■ PU^P 



= tr £ («l|^/ii|n2)(n2|^/i2l"3)"-(n^|^/ijvK)) (2-15) 
niji2,—,nflEA 

which is manifestly gauge invariant. Here, the arbitrary phase factor cancels between \rik) and {rtk\. 
Its gauge invariance is also numerically checked in the lattice QCD Monte Carlo calculation. 
From Tr W P (R, T) on the RxT rectangular loop, we define Dirac-mode projected potential, 

V P (R) = - lim -ln{Tr W P (R,T)}. (2.16) 



On a periodic lattice of V = U x 7V t , we define the Dirac-mode projected Polyakov loop [|1 1J: 
(Lp-) = ±-Tr f\^= i-Tr {U P ) N < = -U £ (m l^kX^fe) • ■ ■ K, |fA>i), (2.17) 

i=l ni,..,njv ( e4 

which is also manifestly gauge-invariant. 



3. Analysis of confinement in terms of Dirac modes in QCD 

In this paper, we mainly consider the removal of low-lying Dirac modes, i.e., the IR-cut case. 
Using the Dirac-mode expansion and projection method, we calculate the IR-Dirac-mode-cut Wil- 
son loop Tr W P (R, T), the IR-cut inter-quark potential V P (R), and the IR-Dirac-mode-cut Polyakov 
loop (Lp}ir in a gauge-invariant manner [JTTJ] . Here, we can directly investigate the relation between 
chiral symmetry breaking and confinement as the area-law behavior of the Wilson loop, since the 
low-lying Dirac modes are responsible to chiral symmetry breaking. 

As a technical difficulty, we have to deal with huge dimensional matrices and their products. 
Actually, the total matrix dimension of (m\U^\n) is (Dirac-mode number) 2 . On the L 4 lattice, the 
Dirac-mode number is L 4 x N c x 4, which can be reduced to be L 4 x N c , using the Kogut-Susskind 
technique [Q]. Even for the projected operator, where the Dirac-mode space is restricted, the matrix 
is generally still huge. At present, we use a small-size lattice in the actual lattice QCD calculation. 
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We use SU(3) lattice QCD with the standard plaquette action at j3 = 5.6 (i.e., a ~ 0.25fm) on 
6 4 at the quenched level. The periodic boundary condition is imposed for the gauge field. We show 
in Fig.2(a) the spectral density p(A) of the QCD Dirac operator D. The chiral property of D leads 
to p(— A) = p(A). Figure 2(b) is the IR-cut Dirac spectral density Pir(A) = p(A)0(|A| — A IR ) with 
the IR-cutoff Air = 0.5a" 1 ~ 0.4GeV for the Dirac eigen-mode. Note that, using the eigenvalue 



X„, the quark condensate (qq)iR with the IR-cut Aa 1r is expressed as [1 1] 



(w)a ir = -^ I J^2- (3-D 

The chiral condensate is largely reduced as (qq)h m / (qq) — 0.02 by removing the low-lying Dirac 
modes in the physical case of m q ~ 5MeV, as shown in Fig.2(c). 




Figure 2: (a) The Dirac spectral density p(A) in lattice QCD at /3=5.6 and 6 4 . The volume V is multiplied. 

(b) The IR-cut Dirac spectral density Pir(A) = p(A)0(|A| - Air) with the IR-cutoff Air = 0.5GT 1 ~0.4GeV. 

(c) The lattice result of (<?<?) a ir / (qq) in the case of IR cut Air = 0.5a -1 , plotted against the current quark 
mass m. A large reduction of (qq)\ lK / (qq) — 0.02 is found in the physical case of m ~ 0.006a _1 ~ 5MeV. 



Figure 3 shows the IR-Dirac-mode-cut Wilson loop (W P (R,T)) = TrW p (R,T), the IR-cut 
inter-quark potential V P (R), and the IR-Dirac-mode-cut Polyakov loop (Lp)jr, after the removal 
of the low-lying Dirac modes. These Dirac-mode projected quantities are obtained in lattice QCD 
with the IR-cut of Pir(A) = p(A)0(|A| - A m ) with the IR-cutoff A K = 0.5a" 1 ~ 0.4GeV. 



4 5 6 
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Re (Lp) m 



Figure 3: The attice QCD results after the removal of low-lying Dirac modes Jl which gives Pir(A) = 
p(A)0(|A| - Air) with the IR-cutoff Ai R = 0.5a- 1 ~ 0.4GeV. (a) The IR-cut Wilson loop Tr W P (R,T) 
(circle) after removing the IR Dirac modes, plotted against RxT. The slope parameter a p is almost the 
same as that of the original Wilson loop (square), (b) The IR-cut inter-quark potential (circle), which is 
almost unchanged from the original one (square), apart from an irrelevant constant, (c) The scatter plot of 
the IR-Dirac-mode-cut Polyakov loop (L/>)ir: its zero-value indicates Z3-unbroken confinement phase. 
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Remarkably, even after removing the coupling to the low-lying Dirac modes, which are re- 
sponsible to chiral symmetry breaking, the IR-Dirac-mode-cut Wilson loop obeys the area law as 
(W P (R,T)) oc e~ aPRT , and the slope a p , i.e., the string tension, is almost unchanged as a p ~ a. 
As shown in Fig.3(b), the IR-cut inter-quark potential V P (R) is almost unchanged from the original 
one, apart from an irrelevant constant. Also from Fig.3(c), we find that the IR-Dirac-mode-cut 
Polyakov loop is almost zero, (Lp)ir ~ 0, which means Z3-unbroken confinement phase. In fact, 
confinement is kept in the absence of low-lying Dirac modes or the essence of chiral symmetry 
breaking [11]. This result seems consistent with Gattringer's formula [Q] and Lang's result [Ojl. 

We also investigate the UV-cut of Dirac modes in lattice QCD, and find that the confining force 



is almost unchanged by the UV-cut [11], as shown in Fig.4. This result seems consistent with the 



lattice result of Synatschke-Wipf-Langfeld [15]. Furthermore, we examine "intermediate(IM)-cut" 



of Dirac modes, and obtain almost the same confining force Ql 1Q, as shown in Fig.5. 




Figure4: (a) The UV-cut Dirac spectral density puv(^) = p(A)0(Auv - |^-|) with Auv = 2a ~ 1.6GeV. 
(b) UV-cut Wilson loop TiW p (R, T) (circle) after removing the UV Dirac modes, plotted against RxT. The 
slope <j p is almost the same as that of the original Wilson loop (square), (c) The UV-cut inter-quark potential 
(circle), which is almost unchanged from the original one (square), apart from an irrelevant constant. 



From these lattice QCD results, there is no specific region of the Dirac modes responsible to 
confinement. In other words, we conjecture that the "seed" of confinement is distributed not only 
in low-lying Dirac modes but also in a wider region of the Dirac-mode space. 

Our lattice QCD results suggest some independence between chiral symmetry breaking and 
color confinement, which may lead to richer phase structure in QCD. For example, the phase tran- 
sition point can be different between deconfinement and chiral restoration in the presence of strong 
electro-magnetic fields, because of their nontrivial effect on chiral symmetry [p]]. 
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